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ABSTRACT 

For toxin/antitoxin (TA) systems, no toxin has been 
identified that functions by cleaving DNA. Here, 
we demonstrate that RaIR and RalA of the cryptic 
prophage rac form a type I TA pair in which the an- 
titoxin RNA is a frans-encoded small RNA with 16 
nucleotides of complementarity to the toxin mRNA. 
We suggest the newly discovered antitoxin gene be 
named ralA for RaI R antitoxin. Toxin RaIR functions 
as a non-specific endonuclease that cleaves methy- 
lated and unmethylated DNA. The RNA chaperone 
Hfq is required for RalA antitoxin activity and appears 
to stabilize RalA. Also, RalR/RalA is beneficial to the 
Escherichia co// host for responding to the antibiotic 
fosfomycin. Hence, our results indicate that cryptic 
prophage genes can be functionally divergent from 
their active phage counterparts after integration into 
the host genome. 



INTRODUCTION 

Toxin/antitoxin (TA) systems are widespread among 
prokaryotes (1). Five different types of TA systems have 
been characterized, depending on the interaction of the TA 
and the nature of the antitoxin (2,3). For type I systems, an 
RNA antitoxin interacts with the toxin transcript and in- 
hibits translation of the toxic protein (4). The toxins and an- 
titoxins of type II systems interact through direct protein- 
protein binding (4). Type III systems rely upon the direct 
interaction of an RNA antitoxin with the toxin protein (5). 
A type IV designation has been proposed for a TA system in 
which the protein antitoxin does not interact with the toxin 
directly but suppresses the toxicity of the toxin by stabilizing 
its target (6), and a type V designation has been proposed 



in which the proteic antitoxin cleaves specifically the mRNA 
of the toxin to prevent the translation of the toxin (7,8). 

In Escherichia coli, most well-studied TA systems belong 
to type II TA systems, in which the labile proteic antitoxin 
binds to the more stable toxin and inhibits its activity. Type 
I and type III systems are less well studied probably due to 
the delay in the identification of small RNAs (sRNAs) (9). 
However, in recent years, bioinformatic searches for puta- 
tive type I TA systems based on tandem copies of the full 
loci and the presence of transmembrane domains has led to 
the identification of multiple copies of potential type I loci 
within new hosts across 774 bacterial genomes (10). 

The role of TA systems in cell physiology, specifically in 
biofilm formation (11,12), per sister cell formation (13,14), 
the general stress response (15,16), phage inhibition (17,18), 
and differential mRNA decay (19,20) is becoming more 
clear. Type II toxins with mRNA endoribonuclease activ- 
ity usually alter gene regulation by cleaving specific mRNAs 
(21). Type I toxin TisB and type V toxin GhoT cause mem- 
brane damage when overproduced and are stress-response 
elements that are actively involved in persistence (7,8,22). 
Type I toxin Hok utilizes post-segregational killing to stabi- 
lize the Rl plasmid (23), but most of the Hok-like toxins in 
the E. coli K-12 chromosome appear inactive (24). Another 
major group of plasmid-based toxins function as gyrase in- 
hibitors, and homologues of the plasmid RK2 ParE/ParD 
TA system in the Vibrio cholerae genome help maintain the 
integrity of the two chromosomes (25). These results in- 
dicate that chromosomal toxins are functionally divergent 
from plasmid loci. 

Prophages or prophage remnants carrying toxic genes 
have been found to harbor TA systems, and five pairs have 
been reported in E. coli K12: RelE/RelB (26), RnlA/RnlB 

(27) , YpjF/YfjZ (28), Ykfl/YafW (28) and CbtA/YeeU 

(28) . Furthermore, two interesting protein toxins have been 
described in cryptic prophage rac in E. coli K 12. KilR is 
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a toxic peptide that inhibits cell division by inhibiting FtsZ 
(29). It has been suggested that RalR alleviates restriction 
modification possibly to protect the bacterial chromosome 
when recombination generates unmodified products by the 
same mechanism as Ral of phage lambda (30,31). How- 
ever, the amino acid identity between the Ral protein in 
phage lambda and the RalR protein in rac prophage is very 
low (24%) (30) and the function of RalR remains unclear. 
Lambdoid prophage rac has lost about 60% of its original 
DNA (32), and ralR is differentially regulated in the devel- 
opment of E. coli biofilms (33). Also, ydaC, downstream of 
ralR, has been named as rcbA for its ability to reduce the 
frequency of double-strand chromosome breaks (34). How- 
ever, it remains unclear whether rcbA encodes RNA or pro- 
tein (34), or whether ydaC is transcribed. 

In this study, unexpectedly, we found that RalR does 
not inhibit restriction modification but instead is part of 
a toxin/antitoxin system and functions as a non-specific 
DNase. Furthermore, we found that the adjacent gene prod- 
uct, a ^raw^-encoded sRNA RalA (for RalR antitoxin), 
functions as an antitoxin for RalR. The ralR and ral A genes 
are adjacent but in the opposite orientation, and RalA RNA 
has 1 6 nucleotides of complementarity to the coding region 
of RalR mRNA. We show that RalA RNA interacts with 
the mRNA of RalR via base-pairing, thus preventing the 
translation of RalR. The activity of antitoxin RalA requires 
RNA chaperone Hfq. Thus, RalR/RalA belongs to a type I 
TA system where the antitoxin sRNA interferes with trans- 
lation of the toxin mRNA via the 16-nt base-pairing. 



MATERIALS AND METHODS 

Bacterial strains, plasmids and growth conditions 

The isogenic E. coli BW25113 K-12 strains and plasmids 
used in this study are listed in Table 1. For deleting and 
overexpressing single genes, the Keio collection (35) and the 
ASKA Hbrary (36) were used. The AralRA mutant strain 
was created using the X Red method (37) using primers 
DralRA-f and DralRA-r (Supplementary Table SI). The 
kanamycin resistance cassette from the Keio and AralRA 
mutants was removed as previously described to ensure 
that only the impact of the deleted gene was studied (38). 
The AralR and AralRA mutations were verified by ampli- 
fying chromosomal DNA using the primers CralRA-f and 
CralRA-v, or CralRA-f2 and CralRA-v2, respectively, and 
sequencing the resulting fragments using CralRA-f (Supple- 
mentary Table SI) to also verify that genes upstream and 
downstream of the mutation (recT and ralR or ydaQ, re- 
spectively) were left intact. Removal of hfq was verified us- 
ing the same method using primer pairs hfq-up and hfq- 
r, hfq-up and hfq-down, hfq-f and hfq-r, or hfq-f and hfq- 
down as described previously (35). All experiments were 
conducted at 37°C in Luria-Bertani medium (LB) (39) un- 
less specified otherwise. Chloramphenicol (30 |xg/ml) was 
used for maintaining pCA24N-based plasmids, and ampi- 
cillin (50 jxg/ml) was used to maintain plasmid pB AD-ralA. 
nullnull 



Cloning of pCA24N-m/i?-m/^, pCA24N-m/i?-m/^- 
G289del and pBAD-m/^ 

The full-length gene region ralR-ralA and partial gene re- 
gion ralR-ralA-G2S9 were cloned as described previously 
(36). They were polymerase chain reaction (PCR) ampHfied 
from E. coli K12 BW25113 chromosomal DNA using one 
front primer (p-ralR-f) and two different rear primers (p- 
ralA-v and p-ralR-ralA-G2S9dQ\-T) (Supplementary Table 
SI). PCR products were phosphorylated using a kination 
kit (Takara, Dalian, China) and purified using a PCR prod- 
uct purification kit (Qiagen, Hilden, Germany). The puri- 
fied products were then ligated into vector pCA24N that 
was digested with the Stul restriction enzyme (New Eng- 
land Biolabs, Beverly, MA, USA). The Hgation mixture was 
transformed into BW251 13 and the Arac mutant. The con- 
structs were verified by DNA sequencing using primer pair 
of p-f and p-r (Supplementary Table SI). Plasmid pBAD- 
ralA was constructed with the full-length gene ral A, and 
124 bp upstream of ral A was controlled by the arabinose 
promoter using primer pair pBAD-ra/v4-f and -r (Supple- 
mentary Table SI) and cloning into A^C6>I-7/z>2dIII-digested 
pBAD/Myc-HisA (Invitrogen, Carlsbad, CA, USA). 



Mutant library creation by epPCR 

Error-prone PCR (epPCR) was conducted on plas- 
mid pCA24N-ra/i^ and pCA24N-ra/i^-ra/y4 using primers 
epPCR-f and epPCR-r (Supplementary Table SI) as de- 
scribed earlier (40). The epPCR program was as follows: 
94°C for 5 min, 30 cycles of 1 min at 94°C, 1 min at 
55°C, and 2 min at 72° C, followed by 10 min at 72°C fi- 
nal extension. The PCR products were gel-purified and di- 
gested using BseRl and HindWl before being ligated into 
pCA24N-ra/i^-ra/v4. The ligation mixture was transformed 
into BW25113 Arac. 



Single and multiple site-directed mutagenesis 

Single site-directed mutagenesis (16) was used to mutate 
the putative start codon of ydaC from ATG to ACG us- 
ing primer pair p-ra/i^-T189C-ra/v4-f and -r, and to intro- 
duce a stop codon TAA into the ninth putative coding re- 
gion of ydaC in pCA24N-ra/i^-ra/^ using primer pair p- 
ralR-ralA-Cl6lA-f and -r (Supplementary Table SI). For 
the 16-nt repeat in RalA, positions 63 (T to C) and 64 (G 
to C) were mutated using corresponding primer pairs. In 
addition, the AGC site in ralA was mutated to AAC us- 
ing p-ralR-ralA-GS3A-{ and -r. Site-directed mutagenesis at 
multiple sites was applied to mutate the GAGC sequence 
into CAAC using p-ra/i^-ra/yl-G81CG83A-f and -r. More- 
over, p-ralR-ralA-Ml-{ and -r were used to mutate the 16- 
nt repeat sequence of ralA on pCA24^-ralR-ralA without 
changing the potential amino acid sequence. p-ralR-Ml-f 
and -r were used to mutate the 16-nt repeat in ralR with no 
amino acid changes in pCA24N-ra/i^ and pCA24N-ra/i^- 
ralA, and p-ra/7^-M2-f and -r were used to mutate the 16- 
nt repeat in ralR in pCA24N-ra/i^ and pCA24N -ral R-ral A 
with amino acid changes (Table 1). The correct mutations 
were verified by DNA sequencing. 
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Table 1. Bacterial strains and plasmids used in this study 



Bacterial strains/Plasmids 


Description 


Source 


E. coli K12 BW25113 strains 






Wild-type 


rmB3 AlacZ4787 hsdR514 A(araBAD)567 A(rhaBAD)568 
rph-1 


(35) 


Arac 


Whole prophage rac removed 


(38) 


ixnjq 


Ahfq Akcin 


This study 


A rnJJ? 


A vrtl /? A Irrtvi 
L\rCllI\ LxhiCiri 


1 Ills aLUUy 


A vnl T? A 


Aval R A Akcin 


This study 








pL.Az4iN 


Cm , Laci^ 


(36) 


pCA24N-ra//< 


Cm^; lacr^, VTs-iac - .^alR 


(36) 


pCA24N-ra/^ 


Cm^; lacl^, Vjs-iac'-'-ydac 


(36) 


pCA24N-ra/i?-raM 


Cm^; lacn, FT5-lac--ralR-ydaC 


This study 


pBAD-ra/^ 


Amp^; FsAD'-'-ralA'^ 


This study 


pCA24N-/z/^ 


Cm^; lacl^, ^TS-lac'-'-hfy 


(36) 


Mutated pCA24N-ralR-ralA plasmids 




This study 


pCA24N-ra/7?-ra/^-G81CG83A 


G to C at position 81 in ralA; G to A at position 83 in ralA 




pCA24N-ra/i?-ra/^-G83A 


G to A at position 83 in ralA 




pCA24N-ra/i?-ra/^-G289del 


Sequences removed after G at position 289 in ralR-ralA 




pCA24N-ra/7?-Tl 89C-ra/^ 


T to C at position 189 in ralR, changed putative ydaC start aa 
from Met to Thr, no change in RalR aa sequence 




pCA24N-ra/i?-ra/^-C161A 


C to A at position 161 in ralA, introduced a stop codon TAA 
at the ninth aa of putative YdaC 




pCA24N-ra/i?-ra/^-G64C 


G to C at position 64 in ralA 




pCA24N-ra/7?-ra/^-T63C 


T to C at position 63, no change in aa sequence 




pCA24N-ra/7?-ra/^-Ml 


Modified the 16-nt repeat in ralA to 

AGCCTCCTTCTTACCT, no change in aa of putative YdaC 




pCA24N-ra/7?-M2-ra/^ 


Modified the 16-nt repeat in ralR, change aa from 

CX\^T Cot- T Arc /^1n Alo (~l\^T Cot- Alo Tlo Alo r^ln 

kjiy-oer-^jiu-J-^ys-o^iu-A-ia to o^iy-t3er-/\ia-iie-/\ia-wiu 




pCA24N-ra/i?-M 1 -ralA-M 1 


Modified the 16-nt repeat both in ralR and ralA to 
AAGTGAGAAGGAGGCT and AGCCTCCTTCTTACCT 




Mutated pCA24N-ra/i? plasmids 




This study 


pCA24N-ra/i?-A154G 


Encodes RalR-mutant RalR K52G 




pCA24N-ra/i?-Ml 


Modified the 16-nt repeat in ralR to 
AAGTGAGAAGGAGGCT, no change in aa 




pCA24N-ra/7?-M2 


Modified the 16-nt repeat in ralR, change aa from 
Gly-Ser-Glu-Lys-Glu-Ala to Gly-Ser-Ala-Ile-Ala-Glu 





Note. Cm^ and Amp^ indicate chloramphenicol and ampicillin resistance, respectively, and aa indicates amino acid. 



Protein expression and purification 

RalR, RalR-mutant, and Hfq each with six histidines and 
10 additional cloning artifacts at the N-terminus were pu- 
rified via BW25113 with pCA24N-ra/i^, pCA24N-ra/i^- 
A154G (encodes inactive RalR K52Q) and pCA24N-/?/^ 
(36). Strains were grown in LB with chloramphenicol and 
were induced with 1 mM IPTG at a turbidity of 0.1 for 
5 h. Cells were collected and resuspended in 25 ml lysis 
buffer [50 mM potassium phosphate buffer (pH 8.0), 300 
mM NaCl, and protease inhibitor cocktail (Sigma- Aldrich, 
USA)]. Then samples were sonicated using a Sonic Dis- 
membrator (Ningbodongzhi, China) at level 2 for 5 min 
twice with 2 s sonication and 4 s break while being cooled 
in ice. Ni-NTA resin (Qiagen) was used according to the 
manufacturer's protocol. Purified proteins were desalted by 
passage on disposable Sephadex G-25 pre-packed PD-10 
columns pre-equilibrated in 20 mM Tris-HCl buffer (pH 
8.0), and the protein concentration was measured by using 
a Bi Yuntian BCA assay kit (Haimen, China). 



DNA cleavage assay 

The DNA cleavage assay was performed as described earlier 
(41,42) with some modifications. Purified RalR and RalR- 



mutant (150 pmol each) were incubated separately with dif- 
ferent substrates at 37°C in 250 mM NaCl, 100 mM Tris- 
HCl (pH 7.5), 10 mM CaCls and 10 mM MgCh for 30 min 
and 120 min. Plasmid pBR322 isolated from BW251 13, ge- 
nomic DNA of BW25113, as well as methylated and un- 
methylated lambda DNA (dam~, dcm~) (1 (xg) were used 
as substrates, respectively. The reaction was stopped by the 
addition of a stop solution (25% glycerol, 0.5% SDS, 0.05% 
bromophenol blue and 50 mM EDTA) and was analyzed 
by electrophoresis on 1% agarose gels stained with SYBR 
safe (Invitrogen). The equal amount of deoxyribonuclease 
I (DNase I; New England Biolabs, Beverly, MA, USA) was 
used as a positive control, and excess EDTA (20 mM) was 
added in the reaction to inhibit the nuclease activity. To 
study whether Ca^^ and Mg^^ are cofactors for RalR, 250 
mM NaCl and 100 mM Tris-HCl (pH 7.5) buffer with or 
without the addition of 10 mM Ca^^ or 10 mM Mg^^ were 
used to conduct the DNA cleavage assay. The RalR-mutant 
was used as a negative control. 

In vitro transcription and RNA labeling 

For in vitro synthesis of T7-RalA sRNA and Tl-ompA 
mRNA, PCR products were obtained from genomic DNA 
using primer pairs T7-ompA-f and -r and T7-ralA-{ and -r. 
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respectively (Supplementary Table SI). The PCR-amplified 
products were used as templates for in vitro transcription 
with T7 RNA polymerase. The T7 RNA polymerase pro- 
moter sequence was included in the forward primers. About 
1 |jLg gel-purified PGR products were used as the templates 
for in vitro RNA reactions with the T7 High Yield RNA 
Synthesis Kit (New England Biolabs, Beverly, MA, USA). 
About 50 pmol of the in vitro transcribed RNA was labeled 
using the Pierce^^ RNA y end biotinylation kit (Thermo 
Scientific, Hudson, NH, USA) according to the manufac- 
turer's protocol. 

RNA cleavage assay 

The RNA cleavage assay was performed as described earlier 
(7). The reaction mixture for the endoribonuclease cleavage 
assay (10 |jl1) contained 2 (xg in vitro synthesized Tl-ompA 
mRNA, 50 mM Tris-HCl (pH 8.5), 100 mM KCl, 2.5 mM 
MgCl2 and 30 [xg of purified RalR or RalR-mutant pro- 
tein. The reaction mixture was incubated at 37°C for 10, 20 
and 30 min and quenched by the addition of an equal vol- 
ume of 2 X TBE-urea sample buffer (Invitrogen). Inactive 
RalR-mutant was used as a negative control. The reaction 
products were resolved by 15% TBE-urea gels. RNA was 
stained with SYBR safe. 

Tricine-SDS-PAGE and western blot analysis 

To investigate RalR protein levels, Tricine-SDS-PAGE and 
a western blot were performed. BW25113 Arac strains 
containing pCA24N, pCA24N-ra/i^, pCA24N-raM and 
pCA24N-ra/7^-ra/yl were grown to a turbidity of 0.2 in LB 
with chloramphenicol, then 1 mM IPTG was added to pro- 
duce RalR for 5 h, and cells were washed with TE buffer. 
Samples were sonicated and total protein was quantified as 
mentioned above. Protein was denatured at 95°C for 5 min. 
The Tricine-SDS-PAGE was performed as described earlier 
(43). Total protein (25 ixg) of each sample was loaded for 
SDS-PAGE, and a western blot was performed using 2.5 
|xg protein of each sample with primary antibodies raised 
against a His-tag (Cell Signahng Technology, Danvers, MA, 
USA) and horseradish-peroxidase-conjugated goat anti- 
mouse secondary antibodies (Bio-Rad, Richmond, CA, 
USA). 

RNA-protein electrophoretic mobility shift assay (EMS A) 

Biotinylated RalA sRNA (4 x 10~^ pmol) was mixed with 
different concentrations of Hfq protein in 20 |jl1 buffer con- 
taining 10 mM HEPES (pH 7.3), 20 mM KCl, I mM 
MgCl2, 1 mM DTT and 5% glycerol at 37°C for 30 min. 
Samples from the reactions were loaded onto a pre-run 6% 
native polyacrylamide gel, and electrophoretic transfer was 
used to bind the proteins to nylon membranes (400 mA for 
30 min) after the separation. After crosslinking for 5 min 
using UV-light, biotin-labeled RNA was detected using a 
LightShift® Chemiluminescent RNA EMSA kit (Thermo 
Scientific, Hudson, NH, USA) according to the manufac- 
turer's protocol. In this assay, in vitro synthesized ompA 
mRNA was used as a negative RNA control and inactive 
Hfq protein was used as a negative protein control. 



Quantitative real-time reverse-transcription PCR (qRT- 
PCR) 

Total RNAs were isolated as described previously (11) and 
was used as the template for the qRT-PCR reaction using 
the Superscript'^^ III Platinum SYBR® Green One-Step 
qRT-PCR Kit (Invitrogen). Primers for qRT-PCR ofralR-, 
ralA- and ra/^-related fragments are listed in Supplemen- 
tary Table SI. The level of rrsG transcript was used to nor- 
malize the gene expression data. 



Primer extension 

5^ end FAM dye (6-carboxyfluorescein) labeled primers 
FAM-ra/^-r2 and FAM-ra/v4-r3 (Supplementary Table SI) 
were ordered from Invitrogen. A total of 30 fxg RNA was 
added to 2 x 10""^ pmol of 5^ end labeled primer, and the 
mixture was added to 3 |xl of lOx first-strand and 37.5 U 
AMV reverse transcriptase (Promega, Madison, WI, USA). 
The RNA mix was annealed to the primers by incubating 
at 37° C for 1 h. To enhance the concentrations of FAM- 
labeled reverse transcriptions, the incubated mixture was re- 
transcribed by adding 6 |jl1 of lOx first-strand buffer and 
75U AMV reverse transcriptase and re-incubated at 37°C 
for another 1 h. The products were screened by an ABI3730 
DNA Analyzer (Applied Biosystems, Foster City, Califor- 
nia, USA) and the results were analyzed using Genemap- 
per (Version 4.1). In this study, only RNA isolated from 
BW25113 Arac/pCA24N-ra/i^-ra/^ induced with 1 mM 
IPTG for 20 min was used due to the low expression level 
ofraMin BW25113. 



Prediction of the valA Rho-independent terminator 

To determine whether RalA sRNA contains a Rho- 
independent structure, we folded its sequence with RNAS- 
TRUCTURE 4.6 (44). A Rho-independent terminator was 
predicted based on properties characterized previously (45). 



Fosfomycin resistance assay 

Metabolic activities in the presence of fosfomycin and glu- 
cose 6-phosphate (to enhance the activity of fosfomycin) 
were measured using reagents from Biolog, Inc. (Hayward, 
CA, USA). Cells were grown to OD600 1-0, then diluted to 
a turbidity of 0.07 in IF-lOa (Cat. No. 72264), and further 
diluted 200-fold into a solution containing IF- 10a, BioLog 
Redox Dye D (Cat. No. 74224), and rich medium (2.0 g/1 
tryptone, 1.0 g/1 yeast extract and 1.0 g/1 NaCl) to a final 
OD600 of 0.0003. Volumes of 100 |jl1 of this cell suspension 
(with 50 mg/ml glucose 6-phosphate and 0 or 0.2 |xg/ml 
fosfomycin) were transferred into 96-well microtiter plates 
and incubated at 37°C. Metabolic activity was measured via 
the absorbance at 590 nm, which indicates the intracellular 
reduced state due to formazan (purple) formed from a tetra- 
zolium dye. Growth in the presence of fosfomycin was deter- 
mined in minimal medium; cells were washed once with M9 
minimal medium with 0.4% glucose (M9-Glu) (46), then di- 
luted to OD600 0.01 in M9-Glu with or without 2 fxg/ml 
fosfomycin. Turbidities were measured at 600 nm. 
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RESULTS 

RalR is a proteic toxin with endonuclease activity 

RASTA-bacteria has predicted RalR of prophage rac to be 
a small toxin (64 aa) (Figure lA, Supplementary Figure SI) 
(47). Indeed, induction of ralR via pCA24N-ra/i^ in E. coli 
results in growth inhibition (Figure IB), along with a 10^- 
fold reduction in colony-forming units (CFU) per milliliter 
in the wild-type host (Figure IC). For the strain with the 
complete rac prophage of 23 kb removed by natural exci- 
sion, Arac, which is devoid of all rac prophage genes (38), 
overproducing RalR was also very toxic (Figure IB and C). 
Hence, RalR is a potent toxin. 

We investigated the cause of the RalR toxicity and found 
it was due to DNase activity. As shown by an in vitro DNA 
cleavage assay, purified RalR cleaves lambda DNA in a pat- 
tern similar to DNase I (Figure 2A). DNase I is an en- 
donuclease that non-specifically cleaves DNA to release di- 
, tri- and oligonucleotides with 5^-phosphorylated and 3^- 
hydroxylated ends (48). After 30 min, ~70% of the lambda 
DNA was cleaved by RalR, and all of the lambda DNA was 
cleaved by DNase I. After 120 min, all of the lambda DNA 
was cleaved by RalR (Figure 2 A and B). Moreover, RalR 
also cleaved native genomic DNA of E. coli K-12 BW251 13 
and plasmid DNA pBR3222 (Supplementary Figure S2A). 

To test whether the endonuclease activity of RalR de- 
pends on modification of substrates by methylation, un- 
methylated lambda DNA was also tested, and the results 
show that RalR cleaves unmethylated DNA equally well, in- 
dicating that methylation is not required (Figure 2B). More- 
over, Ca^^ and/ or Mg^^ were essential for RalR endonucle- 
ase activity (Figure 2C), and as expected, the addition of ex- 
cess EDTA inhibited the activity of both RalR and DNase 
I (Figure 2 A and B). 

Some nucleases degrade all nucleic acids, i.e. both RNA 
and DNA (49). Thus, the endoribonuclease activity of RalR 
was also tested, but no activity was detected in the in vitro 
RNA cleavage assay (Supplementary Figure S2B). There- 
fore, RalR is a non-specific endonuclease which cleaves 
DNA but not RNA. Furthermore, we found that overpro- 
duction of RalR after 5 h leads to filamentous growth as 
a result of the SOS response induced by degrading DNA 
(Supplementary Figure S2C). 

We further determined the key residues in RalR using ep- 
PCR using pCA24N-ra/i^ as the template. After screening 
~400 colonies for toxicity, 39 of them with reduced toxic- 
ity were sequenced, and changes at amino acid position 14, 
37, 39, 49, 51 and 52 reduced RalR toxicity (Supplementary 
Figure S3). One of the mutated proteins (at aa 52) was also 
purified (RalR-mutant, RalR K52Q) (Supplementary Fig- 
ure S4), and since it did not cleave DNA under the same 
conditions (Figure 2 A and B), the DNase activity of RalR 
is not due to contamination. 

DNA associated with j^/aC functions as an antitoxin 

Next, we tested whether previously annotated YdaC works 
as the antitoxin for RalR. YdaC is annotated as a rac- 
encoded protein (69 aa), and the start codon of ydaC over- 
laps with the stop codon of ralR (http://ecogene.org) (Sup- 
plementary Figure SI). The chromosomal region from the 



translational start of ralR to the translational stop of ydaC 
was used to construct a new plasmid pCA24N-ralR-ralA 
(Figure lA). Results showed that when co-expressed with 
ralR, ralA, via the pCA24N-ralR-ralA construct, reduces 
the toxicity of RalR, demonstrating that the DNA associ- 
ated with ydaC is necessary for repression of RalR toxicity 
(Figure IB and C). However, no ydaC-Qncodcd protein was 
detectable via SDS-PAGE or western blot when ydaC was 
tagged with an N-terminal hexahistidine sequence (Figure 
3A). 

To further test whether YdaC protein is the antitoxin, we 
introduced a stop codon by a single nucleotide change into 
ydaC at the putative amino acid position 9 (nucleotide po- 
sition 212 from G to T) (pCA24N-ra/i^-ra/v4-C161A) and 
tested its impact on cell growth. However, the G212T substi- 
tution did not affect the ability of the j<iaC-associated DNA 
to block the toxicity of RalR in cell growth (Figure 3B). In 
addition, when the putative start codon ATG at position 
189 relative to the first nucleotide of start codon of ralR 
(lies inside the coding region of RalR) was mutated into a 
non-start codon ACG (pCA24N-ra/i^-T189C-ra/^), which 
would prohibit the translation of the hypothetical YdaC, the 
DNA associated with ydaC still functioned as an antitoxin 
(Figure 3B). Thus, the DNA associated with ydaC inhibits 
the toxicity of RalR, but ydaC does not encode a protein. 
When we aligned ralR with the region associated with ydaC, 
we found two 16-nt repeats inside ralR and inside ydaC, re- 
spectively (Figure lA). One of the 16-nt repeats lies in the 
coding region of the ralR, and the other, in the opposite ori- 
entation, lies 105-nt downstream of the end of the coding 
region of ralR on the complementary strand. In light of the 
recent discoveries that non-coding regulatory RNAs could 
act as antitoxins by base-pairing with the toxin mRNAs us- 
ing limited complementarity (10-25 nt) (10,50-52), we hy- 
pothesized that the DNA strand complementary to ydaC 
encodes an RNA antitoxin for RalR and propose the name 
RalA for RalR antitoxin. 

RalA is a ^mw^-encoded sRNA 

To further study the antitoxin activity of RalA, we con- 
structed a dual plasmid system to allow expression of ralR 
and ralA independently. Plasmid pCA24N-ra/i^ was in- 
duced by IPTG to express ralR and pBAD-ralA was in- 
duced by L-arabinose to express ralA using the complemen- 
tary strand of the ydaC gQHQ (from 99-nt downstream of the 
stop codon to the start codon of ydaC) (Figure 1 A). pBAD- 
ralA was constructed to have a stop codon (TAA) imme- 
diately after the ATG start codon so that only RNA was 
formed. Toxicity test results indicated that there was a 100- 
fold reduction in viability when ralR is induced by IPTG for 
4 h; however, when L-arabinose was also added to induce 
ralA, the toxicity of RalR was completely inhibited (Figure 
4A). These results show that the RalA sRNA provides the 
antitoxin activity for RalR. 

To map the 5^ end of RalA, we carried out primer exten- 
sion experiments with oligonucleotide FAM-ra/^-r2 with 
complementarity to RalA at regions 34-nt downstream 
from the 16-nt repeat region (Figure 4B). The results re- 
vealed a major extension product that is 1 10 nt in size, which 
indicates that the start of the transcript is 62-nt upstream of 
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Figure 1. RalR is toxic and RalA reduces the toxicity of RalR. (A) The chromosomal region of ralR-ralA in E. coli is shown in the upper panel. ralR is 
shown by the blue arrow, ydaCh shown by the gray arrow, ralA is shown by the red arrow, and the ORFs in the neighborhood of ralR/ralA are shown as 
gray arrows. The direction and the specific region cloned into each of the plasmids used in this study are shown in the lower panel. The numbers indicate 
the position of the related nucleotides (position 1 on the sense strand indicates the first base A of the ralR ORF, and position 1 on the anti-sense strand 
indicates the first base G of RalA sRNA). (B) Cell growth in LB plates supplemented with chloramphenicol (30 |JLg/ml) with and without 0.5 mM IPTG for 
cells producing RalR and RalA in the wild-type host and in the Arac host. (C) CFU test in LB medium supplemented with chloramphenicol (30 |jLg/ml) 
with 1 mM IPTG (added at ODgoo 0-1) in the wild-type host and in the Arac host. Three independent cultures of each strain were evaluated. 



the 16-nt repeat region (Figure 4C, upper panel). No sig- 
nal was detected using another primer that binds further 
downstream (FAM-ra/y4-r3) for primer extension (Figure 
4C, lower panel). 

Since non-coding sRNAs frequently terminate transcrip- 
tion via Rho-independent transcription terminators (53), 
we first searched for a potential Rho-independent termi- 
nation signal in ralA. After folding RalA RNA by RNAS- 
TRUCTURE 4.6 (44), we identified a Rho-independent ter- 
mination signal, with a stem-loop structure composed of 5 
base pairs (GC rich) and a loop of 14 bases followed by 
four stretches of uridine residues (Figure 4D). This stretch 
is likely a Rho-independent terminator since (i) it is pre- 
dicted to form a stem-loop structure, (ii) the formed stem 
has 5-10 base pairs in length and is GC rich with at least 
60% GC base pairs, (iii) the formed loop generally contains 
3-8 bases, (iv) in most structures, the average free energy 
calculated for a Rho-independent terminator is around —7 
kcal/mol and here it was found to be —6.8 kcal/mol, and 
(v) the stem structure usually is followed by a stretch of Us 
(the average number is 4) and here it was found to have 4 
(45). The Rho-independent termination signal thus gives 



an RNA product of 179 nucleotides (Figure 4D), which in- 
dicates that RalA sRNA is wholly contained within what 
was previously annotated as ydaC (on the opposite DNA 
strand). This explains why the ydaC DNA fragment was 
found to reduce RalR toxicity (Figure IB and C). More- 
over, the 5^ end of RalA also formed a stable stem-loop at 
the transcription start site determined by primer extension, 
which helps to stabilize the sRNA molecule (Figure 4D). 
Furthermore, both transcripts were detected in the expo- 
nential phase (OD600 0.8) and stationary phase ( OD600 6.2) 
of the wild-type strain using qRT-PCR with primer pair 
ralR-f and ralR-v (for ralK) and primer pair ralA-f2 and 
ralA-v3 (for ralA) (Supplementary Table SI), which shows 
that both native promoters function in vivo. 

We further confirmed the length of the RalA transcripts 
using qRT-PCR and using cDNA templates synthesized 
from strain- specific primers in exponentially growing cells 
of BW25113 and BW25113 Arac carrying pCA24N-ra/i^- 
ralA (Figure 4B). To test further whether RalA is an un- 
translated RNA, we searched for open reading frames in 
which a putative initiation codon (ATG, TTG or GTG) was 
preceded at an appropriate distance by a ribosome-binding 
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Figure 2. RaIR cleaves DNA. (A) RaIR cleaves lambda DNA after 30 min and 120 min. (B) RaIR cleaves unmethylated lambda DNA (dam~, dcm~) after 
30 min and 120 min. In (A) and (B), the positive control is degradation of DNA by DNase I, and the negative control is the inactive RalR-mutant (RaIR 
K52Q). EDTA blocks RaIR and DNase I activity. (C) RaIR activity requires co-factor Mg^^ and/or Ca^+. 
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Figure 3. RalR is a protein and RalA functions as RNA. (A) Upper panel (SDS-PAGE) and lower panel (western blot) show that there is less RalR protein 
in the presence of RalA. His-tagged RalR (marked with arrows) was produced from pCA24N-ra/i? and pCA24N-ra/i?-ra/y4 and has six histidine residues 
and 10 other amino acids attached at the N-terminus (expected size of ~9.05 kDa). Production of RalR and RalA was induced via 0.2 mM IPTG at a 
turbidity of 0.1 for 5 h, and no IPTG treatment was used for the negative control. (B) Cell growth on LB plates with 0.5 mM IPTG in the Arac host. 
T189C refers to the y<iaCmutation in which the putative start codon ATG at position 189 was mutated into ACG, and C161A refers to the ralA mutation 
containing a stop codon instead of the ninth codon of putative protein YdaC. 



site. This transcript does not contain a potential initiation 
codon preceded by a potential ribosome-binding site, thus 
it is unlikely that RalA encodes a protein. 

RalA functions through base-pairing with RalR 

We next investigated whether RalA RNA functions through 
a base-pairing mechanism. epPCR was employed to screen 
for the mutations in ralR and ralA that result in reduced an- 
titoxin activity using pCA24N-ralR-ralA as template. After 
screening ~400 colonies for toxicity, 20 of them with toxic- 
ity were sequenced, and the results showed that half of them 
had mutations within the two 1 6-nt repeats (Figure 5 A) in 
the ralR gene and in the ralA gene. These results suggest that 
the 1 6-nt repeat is important for RalA RNA to inhibit RalR 
mRNA. 

To further investigate the involvement of the 16 nt in 
the base-pairing of the toxin and antitoxin RNAs, we per- 
formed site-directed mutagenesis on the 16-nt repeats as 
well as other neighboring regions. Our results showed that 
single mutations in the 16-nt repeat of ralA reduced its an- 
titoxin function (pCA24N-ra/i^-raM-G64C and pCA24N- 
ralR-ralA-T 63C), and four mutations in the 16-nt repeat 
of ralA completely abolished the antitoxin activity of RalA 
(pCA24N-ra/i^-raM-Ml). In contrast, a single point muta- 
tion (pCA24N-ra/i^-ra/v4-G83A) and two point mutations 
(pCA24N-ra/i^-ra/^-G81CG83A) that changed one or two 
nucleotides in regions outside of the 16-nt repeat of ralA did 
not affect the antitoxin activity of the RalA RNA (Figure 
5B). In addition, to determine the minimal length of RalA 
RNA required for maintaining the function as an antitoxin, 
we removed all of the sequence after position 289, which 



eliminated the 16-nt repeat in RalA (pCA24N-ra/i^-ra/v4- 
G289del). As a result, the antitoxin function of RalA was 
completely abolished (Figure 5B). These results confirm the 
importance of the 16-nt repeat in the RalR/RalA TA pair 
in regard to the interaction of the toxin and antitoxin. 

Moreover, when the 16 nt in ralR was altered to 
cause a synonymous mutation and non-synonymous mu- 
tation using pCA24N-ra/i^, the synonymous mutation re- 
tained the RalR toxicity (pCA24N-ra/i^-Ml) while the non- 
synonymous mutation led to a RalR variant (pCA24A^- 
ralR-M2, pCA24N-ralR-M2-ralA) that was no longer toxic 
(Figure 5A). These results confirm that RalR functions as a 
protein and, moreover, that amino acid changes at the 16-nt 
repeat coding region affect RalR toxicity. 

To further test whether changes in the ralA 16-nt repeat 
can neutralize the toxic effect of RalR when the correspond- 
ing mutation is encoded by ralR, we introduced two com- 
plementary mutations and constructed pCA24N-ra/i^-Ml- 
ralA-Ml, which contains two identically mutated 16-nt re- 
peats. The results show that the mutated ralA retains its an- 
titoxin function when the mutation matches that of ralR, in 
contrast to the result that mutated ralA could not attenu- 
ate the toxic effect of the unmutated ralR (pCA24N-ralR- 
ralA-Ml) (Figure 5 A). Therefore, we provide several lines 
of evidence that RalA is a novel RNA antitoxin in which 
its activity to block the toxicity of RalR relies on a 16-nt 
repeat. 

RalA likely inhibits RalR mRNA translation 

To determine how RalA RNA attenuates the toxic effects 
of RalR, we checked the RalR protein levels with and with- 
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Figure 4. RalA is a trans-Qncoded sRNA. (A) The dual plasmid system pCA24N-ra/7? and pBAD-ra/^ was used to produce RalR protein and RalA sRNA 
by 1 mM IPTG and 0.1% L-arabinose in the Arac host, respectively The two inducers were added simultaneously at OD600 0-1 and the results of 4 h 
induction with or without the inducers are shown here. Two independent cultures of each strain were evaluated. (B) The position and direction of the 
primers used for primer extension and qRT-PCR to identify the start of transcription of ral A. The number 1 indicates the first base G of RalA sRNA. The 
length of the four fragments (Nos. 1-4) used for qRT-PCR are indicated. The Ct values of each fragment detected in the wild-type strain (ODeoo 1-0) or in 
the Arac/pCA24N-ra/i?-ra/^ strain are indicated after 1 mM IPTG for 20 min once the OD600 reached ~1.0, and the lower Ct indicates higher expression 
level. N/A indicates no signal was detected. (C) Results of the primer extension using FAM-ra/^-r2 or FAM-ra/^-r3, respectively. (D) RNA structure of 
RalA predicted by RNASTRUCTURE 4.6, with a stable start stem-loop at 5' end and a Rho-independent terminator at 3' end. 



out RalA. Through SDS-PAGE and western blot analysis, 
we found that the protein levels of RalR are reduced in the 
presence of RalA (Figure 3A), which indicates that RalA 
interferes with the accumulation of RalR protein. Base- 
pairing between sRNA and its target mRNA usually leads 
to repression of protein levels through translational inhi- 
bition and/or mRNA degradation (51). Our qRT-PCR re- 
sults showed that in the presence of more RalA, the levels 
of RalR mRNA were not reduced substantially (Table 2), 
thus it is unlikely that RalA RNA promotes RalR mRNA 
degradation. Translational inhibition of sRNA with a target 
gene often occurs when sRNA base-pairs with the 5^-UTR 
of target mRNA; however, other locations for base-pairing 
and consequent mechanisms of regulation are possible (51). 

Hfq is required for RalA activity 

To facilitate limited base-pairing for the fra^^-encoded sR- 
NAs and the target mRNAs, in many cases the RNA chap- 
erone Hfq is required for rra^^-encoded sRNA-mediated 



regulation (51). The antitoxin RalA was tested in an hfq 
deletion mutant strain and, as expected, the hfq mutation 
aboHshed antitoxin activity of RalA (Figure 6 A and B). 
Hence, Hfq is required for RalA RNA activity. 

To provide additional evidence for the requirement of 
Hfq for RalA RNA activity, we investigated direct binding 
of Hfq to RalA RNA in vitro. Hfq interacts strongly with 
small non-coding RNAs at single-stranded AU-rich regions 
(54), and RalA is AU rich (63.1%) and has four stretches 
of more than seven continuous U/As. RNA-protein bind- 
ing was performed, and Hfq was found to bind RalA RNA 
(Figure 6C). In contrast, Hfq did not bind to the negative 
control (3^ end of the coding region of OmpA mRNA, Fig- 
ure 6C). Moreover, it seems that Hfq helps to stabilize the 
RalA transcript as RalA RNA was degraded faster in the 
absence of Hfq as shown by qRT-PCR, while both RalR 
and RalA RNAs appeared stable in the wild-type cells with 
half-times of more than 20 min (Table 3). 
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Figure 5. RalA interacts with RalR mRNA through 1 6-nt repeats. (A) The 1 6-nt repeat sequences of both RalR and RalA are shown in the upper schematic, 
and the numbers indicate the positions of the related bases. The table indicates the 16 nt in each plasmid, and the mutated nucleotides are marked in red 
font. 'Toxicity' indicates whether the toxin was toxic under the conditions indicated. Lower figures indicate growth on LB plates with cells producing both 
RalR and RalA via 0.5 mM IPTG in the Arac host. (B) Cell growth on LB plates producing both RalR and RalA via 0.5 mM IPTG in the Arac host. Two 
independent cultures of each strain were evaluated, and no IPTG treatment was used as the negative control. 



RalR increases resistance to fosfomycin 

Large-scale phenotypic screening of all Keio mutants from 
PortEco showed that the l^ralR knockout strain has a neg- 
ative fitness score when exposed to fosfomycin (55). Fos- 
fomycin is a broad-spectrum antibacterial agent that tar- 
gets bacterial mucopeptide synthesis by inhibiting phos- 
phoenolpyruvate transferase, the first enzyme involved in 
the synthesis of peptidoglycan (56). Our results show that 
both the /^ralR and the AralRA strain are more sensitive 
to fosfomycin (0.2 (xg/ml) compared to the wild-type strain 
based on their reduced metabolic activity (Figure 7A and 
B). Corroborating these results, the AralRA strain is more 



sensitive to fosfomycin (2 |xg/ml) compared to the wild-type 
in terms of growth as well (Figure 7C and D). Hence, the 
RalR/RalA TA locus has an impact on cell physiology and 
is beneficial for resistance to an antibiotic. 

DISCUSSION 

Collectively, our results support that RalR/RalA forms a 
type I TA pair. These results are: (i) both genes are small 
and are adjacent, (ii) RalR functions as a toxin that in- 
hibits growth, (iii) RalA functions as an antitoxin and is a 
non-coding RNA, (iv) RalA blocks RalR-mediated toxic- 
ity by the base-pairing of 16 key nucleotides and prevents 
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Table 2. Levels of RalR and RalA transcripts (Ct values ± S.D.) quantified by qRT-PCR in BW251 13/pBAD-ra/^ 
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Figure 6. Hfq is required for RalA antitoxin activity. (A) Cell growth on LB plates supplemented with chloramphenicol (30 |JLg/ml) with and without 0.5 
mM IPTG for cells producing RalR and RalR/RalA in the wild-type host and in the Ahfq host, respectively. (B) CFU test over time for cells producing 
RalR and RalR/RalA in the wild-type host and in the Ahfq host, respectively. Overnight cultures were diluted to ODgoo 0-1, and 1 mM IPTG was added 
initially. This assay was performed twice with two independent cultures, and one standard deviation is shown. (C) EMSA shows that Hfq binds to RalA 
sRNA (left panel) but not to the 3' end of the coding region of ompA mRNA (middle panel). The binding of Hfq to RalA increases with increasing Hfq 
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the translation of RalR mRNA, and (v) Hfq is required 
for RalA antitoxin activity. As a novel type I TA system, 
RalR is the first toxin that functions as an endonuclease in 
E. coli. A recent study showed that when compared with a 
wild-type strain, a ydaC deletion mutant has an increased 
level of double-strand breaks that is remedied by expres- 
sion of ydaC from a plasmid, suggesting that ydaC helps 
to maintain the integrity of the bacterial chromosome (34). 
These results fit well with our results here in that deletion 
of ralA (on the opposite strand as ydaC) leads to activation 
of toxin RalR that should result in double- stranded DNA 



breaks. Although overproduction of RalR leads to filamen- 
tous growth due to RalR-mediated DNA damage that ac- 
tivates the SOS response, it is not clear whether RalR can 
trigger cell death under physiological conditions. 

Another novel feature of this TA system is that RalA is 
a fra^^-encoded sRNA antitoxin that neutralizes the pro- 
tein toxin by base-pairing within the toxin mRNA coding 
region. Traditional type I toxin and antitoxin genes are usu- 
ally located at the same locus, and sRNA antitoxins can 
have greater than 60-nt base-pairing with the 5^ UTR (e.g. 
Ldr of the Ldr/Rdl family in E. coli (57)) or the y UTR 
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Table 3. Levels of RalR and RalA transcripts (Ct values ± S.D.) quantified by qRT-PCR in BW251 13 and BW251 13 ^hfq after the addition of rifampicin 
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Figure 7. The RalR/RalA TA system increases resistance to fosfomycin. (A) Metabolic activity of the wild-type, AralR and the AralRA strains subjected 
to 0.2 |JLg/ml fosfomycin with 50 |JLg/ml glucose-6-phosphate to enhance fosfomycin activity and (B) metabolic activity with 50 |JLg/ml glucose-6-phosphate 
alone (negative control). This assay was performed twice with two independent cultures, and one standard deviation is shown. Growth of the wild-type 
and AralRA strains (C) with 2 |JLg/ml fosfomycin and (D) without fosfomycin. Two independent cultures were evaluated, and one standard deviation is 
shown. 
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of the toxin mRNA (e.g. TxpA/RatA in Bacillus subtilis 
(58)), and the two transcripts can also overlap (e.g. Ibs/Sib 
in E. coli) (50,52). Recently, four unconventional type I TA 
systems have been reported where the antitoxins are tran- 
scribed divergently from the toxin gene and have limited 
complementation regions to the toxin RNA (often 18-21 
nucleotides), including TisB/IstR-1 (59), ShoB/OhsC (50) 
and ZorO/OrzO loci (50). Here with RalR/RalA, we pro- 
vided evidence that RalA reduces RalR toxicity via the com- 
plementarity of 16 nt. For the ^ra/i^-encoded base-pairing 
sRNAs, the potential base-pairing with target mRNA typi- 
cally is 10-25 nt. Studies on ZorO/OrzO showed that for the 
18-nt complementarity region, not all of the base-pairing 
interaction is necessary, and it is dependent on whether the 
complementary region is contiguous or not (52). Our results 
confirmed this and showed that 16 nt of complementarity 
(with 1 1 nt of contiguous perfect complement) is sufficient 
for base-pairing interaction between toxin and antitoxin. 

Cryptic prophage rac is a DNA fossil since it was ac- 
quired over 4.5 million years ago (60), yet it has become 
a stable resident in the E. coli chromosome. Ral in active 
phage lambda and RalR in prophage rac share a low se- 
quence identity of 24%, indicating that the two proteins 
are not related. In contrast, a high identity of 77% ex- 
ists for P22 Ral and lambda Ral with highly conserved N- 
and C-termini (Supplementary Figure S5). Moreover, an 
even higher identity of 94% is found for recombinant En- 
terobacteria phage phi21 Ral and lambda Ral with highly 
conserved N- and less conserved C-termini (Supplemen- 
tary Figure S5). Therefore, Ral is more conserved in phages 
while RalR is more conserved in E. coli strains. Ral of 
lambda phage acts by modulating the restriction and mod- 
ification (RM) activities of the type I restriction systems in 
E. coli (61), and RM systems and TA systems are related 
(3). In general, they originate from phages and plasmids 
and are used to invade hosts or to keep them in the host by 
post-segregational killing. However, after integration into 
the chromosome, the function of these two systems has di- 
verged. We found that RalR in E. coli functions differently, 
acting as a toxin by cleaving DNA, and that it belongs to a 
type I toxin-antitoxin system. We show here that removal 
of ralRA reduces the resistance of the cell to the antibiotic 
fosfomycin; hence, there is a clear benefit to the cell for har- 
boring the RalR/RalA TA system. Therefore, our results 
provide further evidence that the genes of cryptic prophage 
impact cell physiology and that these genes may be used to 
increase resistance to stress (38). 

Besides RalR/RalA, five other TA systems located in the 
E. coli K12 cryptic prophages have been reported, includ- 
ing the well-studied RelE/RelB (in Qin prophage) (26) and 
four recently identified TA pairs YpjF/YfjZ (in CP4-57) 
(28), RnlA/RnlB (in CP4-57) (27), Ykfl/YafW (in CP4-6) 
(28) and CbtA/YeeU (in CP4-44) (6). These TA systems 
are located in prophages that are not inducible by SOS re- 
sponses (38), suggesting a relatively stable residence in the 
host genome. Type I TA pair TxpA/RatA in B. subtilis is 
also in a phage-like element, and the phage-like element is 
retained in the host chromosome for spore formation (58). 
TA systems have been linked to persistence, biofilm forma- 
tion and the general stress responses, and ralR was found 
to be highly induced in the later stage of E. coli biofilm de- 



velopment (24 h biofilms compared to 4 h biofilms) (33). 
Indubitably, other uncharacterized TA systems in cryptic 
prophages might also contribute to the stress response of 
the host, and there are probably additional intricate regula- 
tion mechanisms that need to be investigated. 

Unconventional type I TA systems are more limited in 
distribution than typical type I TA systems, and are found 
mainly in Escherichia, Shigella and Salmonella species. By 
comparing the phylogenetic tree of those identified type I 
loci with the host taxonomy, earlier studies suggest that 
type I loci have not been freely disseminated by horizon- 
tal gene transfer but instead may have a common ancient 
ancestor (10). This agrees with what we found for TA pair 
RalR/RalA in rac prophage. The discovery of a high occur- 
rence of TA systems in E. coli prophages will help to search 
for new TA systems in other prokaryotes. Since TA genes 
from bacterial prophages are actively expressed in the hu- 
man gastrointestinal tract, TA systems may play a dynamic 
role in more ecosystems than expected (62). The mechanism 
of harboring TA systems in prophages, the origins of these 
TA systems, and the regulation of these loci in response to 
various stress conditions need to be investigated further. 
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